from scipy import *
from numpy import *

n = 3
ind1 = arange(n)
t0 = cos(pi*(ind1+1)/n)  
x = arange(n)
def f(p): return (1.0 if p==1 else 0.0)
#print f(x)

xx, yy = meshgrid(x, x)
#print kron(2,2)
def krd(x, y):
  return int(x == y)

vector_krd = vectorize(krd)
a = arange(n)
b = arange(n)
print vector_krd(xx, yy)
#print xx
#print yy
#print f(yy,xx)
def f(l): 
    if l==0: return 1
    else: return 2
#print f(0), f(1), f(2)